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We present an improved laboratory test of Lorentz invariance in electrodynamics by testing the 
isotropy of the speed of light. Our measurement compares the resonance frequencies of two or- 
thogonal optical resonators that are implemented in a single block of fused silica and are rotated 
continuously on a precision air bearing turntable. An analysis of data recorded over the course of 
one year sets a limit on an anisotropy of the speed of light of Ac/c ~ 1 x 10~^^. This constitutes 
the most accurate laboratory test of the isotropy of c to date and allows to constrain parameters of 
a Lorentz violating extension of the standard model of particle physics down to a level of lO"^'^. 



The theory of special relativity formulated in 1905 [2] 
revealed Lorentz invariance as the universal symmetry of 
local space-time, rather than a symmetry of Maxwell's 
equations in electrodynamics alone. This striking insight 
was drawn from two postulates: (i) the speed of light 
in vacuum is the same for all observers independent of 
their state of motion, and (ii) the laws of physics are the 
same in any inertial reference frame. Today, local Lorentz 
invariance constitutes an integral part of the standard 
model of particle physics, as well as the standard theory 
of gravity, general relativity. Still, there have been claims 
that a violation of Lorentz invariance might arise within 
a yet to be formulated theory of quantum gravity [3H5j. 
Given a lack of quantitative predictions, the hope is to 
reveal a tiny signature of such a violation by pushing test 
experiments for Lorentz invariance across the board. An 
overview of recent such experiments can be found in (Hj. 

Previous measurements testing the isotropy of the 
speed of light, often referred to as modern Michelson- 
Morley experiments [TU], have compared the resonance 
frequencies of optical [TTHTl] or microwave [1] [T5] cav- 
ities, which were either actively rotated on a turntable 
or relied solely on Earth's rotation. The most precise of 
these have tested the isotropy of c at an accuracy of a 
few parts in 10^^ limited by relative resonator frequency 
stability. 

The experiment presented here improves on this by one 
order of magnitude, based on an optimized cavity design 
and rotation on a precision turntable that allows to mini- 
mize systematic effects. The basic principle is depicted in 
Figure [T] At the core of the experiment are two crossed 
optical Fabry-Perot resonators. We compare their res- 
onance frequencies by stabilizing two Nd:YAG lasers to 
these cavities and taking a beat note measurement. The 
resonance frequency ^ of a linear Fabry-Perot cavity de- 
pends on the speed of light c along its optical axes as 
given by 

v = mc/2L (1) 

where m is an integer number and L is the length of the 
resonator. Thus, to detect an anisotropy of the speed of 
light Ac ~ Cx — Cy we continuously rotate the setup and 
look for a modulation of the beat frequency lS.v. Since 
the light in the cavities travels in both directions and 



c refers to the two-way speed of light, such an isotropy 
violation indicating modulation would occur at twice the 
rotation rate. 




FIG. 1: Left: High-finesse fused silica resonators used in this 
experiment. Right: Basic principle of the experiment. The 
frequencies of two lasers, each stabilized to one of two or- 
thogonal cavities, are compared during active rotation of the 
setup. (Photograph by E. Fesseler) 



I. THE EXPERIMENT 

The experiment applies a pair of crossed optical high- 
finesse resonators implemented in a single block of fused 
silica (Figurejlj). This spacer block is a 55 mm x 55 mm x 
35 mm cuboid with centered perpendicular bore holes of 
10 mm diameter along each axis. Four fused silica mir- 
ror substrates coated with a high-reflectivity dielectric 
coating at A = 1064 nm are optically contacted to ei- 
ther side. The length of these two crossed optical res- 
onators is matched to better than 2 nm. The finesse of 
each resonator (TEMoo mode) is 380 000, resulting in a 
linewidth of 7 kHz. Two Nd:YAG lasers at A = 1064 nm 
are stabilized to these resonators using a modified Pound- 
Drever-Hall method IF]. Tuning and modulation of the 
laser frequency is achieved with piezo electric actuators 
attached to the laser crystal. Mechanical resonances of 
the piezo-electric actuators at fm = 444 kHz and 687 kHz 
respectively are used for modulation of the laser frequen- 
cies. The light reflected from the cavities is detected and 
demodulated at 3/m to generate an error signal. Ther- 
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FIG. 2; Allan deviation calculated from a comparison of the 
two stabilized laser frequencies (a) while the setup is rotating 
and (b) with a stationary non rotating setup. 



mal effects from dissipation of laser power inside the res- 
onators are minimized by coupling less than 50 /iW opti- 
cal power into the cavities. 

The cavities are installed in a vacuum chamber to en- 
sure light propagation in vacuum inside the resonators 
and to reduce the influence of environmental noise, e.g. 
thermal fluctuations and vibrations in the laboratory. By 
placing the coupling optics inside the vacuum chamber as 
well, a high pointing stability of the laser beams incident 
on the resonators could be achieved. The custom-made 
vacuum chamber features several stages of thermal insu- 
lation and is placed upon an active vibration isolation 
system [EWL, 350-M) (Figure [s]). 

To compare the stabilized laser frequencies, fractions of 
1 mW of each laser beam are split off and are overlapped 
on a fast photodiode to generate a beat note at the dif- 
ference frequency — v\ — vi. By choosing appropriate 
longitudinal modes of the cavities, this frequency is set 
to < 2 GHz and counted with a sampling time interval 
of one second. Since the length of the two cavities is 
defined by a single monolithic block, drifts due to ther- 
mal expansion are largely the same for both resonators 
and thus cancel in a measurement of the difference fre- 
quency. We have observed a reduction of the relative 
drift to below 0.01 Hz/s as compared to lOOHz/s abso- 
lute frequency drift of the individual resonators. 

The stability of the frequency beat note is then char- 
acterized by calculating the Allan deviation. There is 
a pronounced flicker floor from 1 s to 200 s at a level of 
~ 1.5 X 10~^^ (Figure [2| presumably caused by thermal 
noise of the mirror substrates. This agrees with an esti- 
mate of the thermal noise level of our cavity based on a 
model of Numata et al. [T7j . 

To enable continuous rotation of the resonators we em- 
ploy an air bearing turntable which carries the complete 
laser stabilization setup (Figure |3]) and rotates at a cho- 
sen rate of Ttt = 45 s. As opposed to using Earths ro- 



tation alone, such active rotation allows us to perform 
hundreds of rotations per day, while taking advantage of 
the excellent mid-term frequency stability of the cavities 
(Figure|2| . On the other hand, active rotation potentially 
causes a systematic modulation of the beat frequency 
and might thus mimic an anisotropy signal. For exam- 
ple, gravitational or centrifugal forces that act on the 
resonators may get modulated with the turntable rota- 
tion and therefore modulate the length of the resonators. 
However, most of these effects lead to a modulation at a 
rate of Wtt = 27r/Ttt so that they are in principle distin- 
guishable from the anisotropy signal searched for at 2ujtt ■ 
Moreover, if the data spans more than one day, system- 
atic effects with a fixed phase in the laboratory average 
out in the analysis for an anisotropy of c that is fixed rela- 
tive to a sidereal frame. Although such an analysis helps 
to discriminate a sidereal anisotropy signal from system- 
atics, a large effort was still made to reduce systematic 
effects both at 2cjtt and wtt- 

First of all, we use a high precision air bearing 
turntable specified for < 1/irad rotation axis wobble. 
Furthermore, we also prevent long-term variations of the 
rotation axis tilt, caused for example by daily fluctua- 
tions of the building tilt of several /xrad. For this we apply 
an active stabilization [O] that keeps the rotation axis 
vertical to better than 1 /irad, which reduces the effect 
from a periodic deformation of the cavities to frequency 
variations of less than 0.1 Hz in amplitude. Effects from 
varying centrifugal forces are also reduced below an am- 
plitude of 0.1 Hz by an active stabilization of the rotation 
rate. Further measures include balancing the center of 
mass of the table (< 1mm offset from the rotation axis) 
and shielding the lasers and optics outside the vacuum 
chamber against air currents and temperature gradients 
in the laboratory. 

At the chosen rotation period of 45 s these measures 
reduce residual systematic frequency variations at 2a;tt 
to amplitudes below 0.1 Hz. This corresponds to a frac- 
tional frequency shift of Aiy/i^o = 3 x 10~^®, which is well 
below the relative frequency stability of the beat note on 
the timescale of a single rotation (see Figure [2]) . While 
even faster rotation would have allowed to acquire more 
data and thus improve statistics, it resulted in increased 
residual systematic effects presumably due to modulated 
centrifugal forces and was thus not implemented. 

Measurements with this setup have been performed in- 
termittently during a time span of more than one year 
from May 2007 to June 2008. The total data includes 
recordings of the beat frequency, time and rotation an- 
gle at a sampling interval of 1 s from more than 130 000 
turntable rotations. 



II. ANALYSIS FOR AN ANISOTROPY SIGNAL 

In what follows we first give a phenomcnological, i.e. 
largely model-independent, description of an anisotropy 
signal and present results from a corresponding analy- 
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sis. In Section III Al and IIIBI we then use these results 
to determine parameters of two different test theories for 
Lorentz violation. Throughout this analysis, we adopt 
the inertial Sun centered celestial equatorial coordinate 
frame (SCCEF) as used for the analysis of similar, pre- 
vious such experiments |18 . This coordinate system has 
the Z axis pointing north, the X axis pointing in the di- 
rection of the vernal equinox point, and the Y axis such 
that (X,Z,Y) form a right handed set. 

Let us first consider the special case of the apparatus 
located at the North pole with the turntable rotation axis 
aligned with Earth's rotation axis. An anisotropy in the 
equatorial XY-plane {cx 7^ cy) then causes a modula- 
tion of the beat frequency with the rotation of the setup. 
As noted above, this modulation would be at twice the 
rotation frequency, i.e. 2wrot- If we fix the time axis rela- 
tive to some arbitrary instant t = 0, we can describe this 
signal as 



1^1 - V2 



1^0 



= S' sin 2uji-ott + C cos 2ujrott, (2) 



where Pq ~ 282 THz and S",C" - £2c^. The rotation 
of the setup within the sidereal frame of reference is a su- 
perposition of turntable rotation ojtt and Earth's sidereal 
rotation at cu^ such that Wiot — "^e ^ ^tt > pl^^ '^^ minus 
depending on the sense of turntable rotation. Since in 
our experiment we have wtt ^ w®, we can describe the 
anisotropy signal as a fast modulation at 2a;tt 



= S sin 2cjtti + C cos 2wtt^, 



(3) 



with amplitudes S and C that slowly vary with Earth's 
rotation as given by 
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FIG. 3: Schematic of the complete rotating setup (top) and 
the custom-made vacuum chamber (bottom). TS=Tilt Sen- 
sor, PDH=Pound-Drever-Hall laser stabilization electronics. 



C = S" sin 2w®i + C cos 2a;©t. (5) 

This daily modulation is essential to discriminate an 
anisotropy signal from constant or slowly varying sys- 
tematic effects caused by active rotation as described in 
Section |lj Only systematic effects subjected themselves 
to a 23.93 h modulation would mimic such a sidereal 
anisotropy signal. 

Next, we consider an experiment located at an ar- 
bitrary geographical latitude x such that Earth's axis 
and the turntable rotation axis do not coincide anymore. 
While this reduces sensitivity to an anisotropy in the 
equatorial XY-plane, it additionally provides sensitivity 
to an anisotropy in the XZ and YZ-plane. Furthermore, a 
modulation at w© in addition to that at 2a;0 will appear. 
We thus generalize the above expressions to the following 



anisotropy signal (see |18| for a formal derivation) 

S* = S'o + Ssi sin(w(g(t - io)) + S^i cos(w0(t - to)) 
+ 5s2 sin(2a;©(t - to)) + Sc2 cos(2w©(i - to)), (6) 

C = Co + Csi sin(a;©(t - to)) + Cd cos(a;©(t - to)) 
+ Cs2 sin(2a;e(< - to)) + Cc2 cos(2we(^ " ^o)) , (7) 

where the phase is fixed by to chosen in accordance to the 
adopted reference frame conventions. Again, this daily 
modulation at w© and 2uj^ makes it possible to distin- 
guish between a sidereal anisotropy of c and systematic 
effects due to the active rotation. 

To analyze our data for a Lorentz violation signal of 
the above form we proceed in two steps. First we deter- 
mine the modulation amplitudes C and S as modeled in 
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Eq.(j3]) from short samples of the data set spanning 10 ta- 
ble rotations each. These samples extend over 450 s each 
such that we may neglect a possible modulation due to 
Earth's rotation within each sample. We use a fit func- 
tion of the form 

l^vlv = S'sin(2a;tt(i-io)) + C'cos(2a;tt(t-io)) 

-h As sin(wtt(i - to)) + cos(wtt(t - *□)) 
+ ^o + ^ii, (8) 

where Aq and A\ account for an arbitrary offset and a 
linear drift while Ag and Ac account for residual sys- 
tematics at Wtt- In accordance with the reference frame 
convention of [18 the starting time t'^ is determined by 
the first instant of the measurement at which one of the 
two resonators is oriented along the North-South direc- 
tion. We choose a sample size of 10 rotations for each 
fit rather than fitting single rotations to reduce the cor- 
relation of a small linear drift and a sinusoidal variation 
of the beat frequency. We found, however, that choosing 
different sample sizes of n = 2 to 20 does not significantly 
change the final results. 

From each sidereal day (23.93 h) of measurement we 
obtain a distribution of 192 values of S and C, and each 
value is assigned the mean time of the respective data 
sample. In total we obtain a distribution of 13384 values 
for S and C as shown in Figure |4^. 

Next, each 23.93 h interval of these distributions is fit- 
ted with equations ([6]) and ([7| to determine whether there 
is any daily modulation as a consequence of a sidereal 
anisotropy of c. The results are shown in Figure |4]3-|-c as 
well as in Figure [Sj Each graph corresponds to a pair of 
sidereal modulation amplitudes and shows a distribution 
of 64 data points. Each point is determined from one day 
of the measurement. The standard error associated with 
each one-day data point is on the order of ~ 5 x 10~^^. 
If for each sidereal modulation amplitude we take the 
mean value of the corresponding 64 data points, we find 
a standard error on the order of ^ 7 x 10~^^ for each dis- 
tribution and no deviation from zero by more than three 
standard errors. 

Figure [4ji shows the results for the modulation ampli- 
tudes 5*0 and Cq. These are solely connected to a modu- 
lation at 2wtt (Eq.([3|, ([6]) and ([7|) and thus are strongly 
effected by any residual systematic effects fixed to the 
laboratory frame. Single points of these amplitudes de- 
viate by several standard errors from zero, however, over 
the complete measurement span of one year, the data 
points vary in magnitude and phase and thus average 
out. The variation of these amplitudes over time can 
also be seen from the bottom graphs in Figure [5] 

Overall, we conclude that no significant evidence for an 
anisotropy of c fixed relative to a sidereal frame can be 
claimed from this data. This of course assumes otherwise 
uncorrelated noise, e.g. no annual phase shift of a non- 
zero anisotropy signal. 



A. Analysis in the framework of the minimal 
standard model extension 



The above phenomenological results can be further 
evaluated as a test of Lorentz invariance in electro- 
dynamics, adopting the Lorentz violating extension of 
the standard model of particle physics by D. Colla- 
day and V.A. Kostelecky et al. [191 EO]- In this stan- 
dard model extension (SME), Lorentz violation in elec- 
trodynamics is modeled by extending the Lagrangian 
of the photonic sector £ — — ^-F'^^Fp^ with a term 
/:oxt = -\i)iF)t.^n\F^'''F>^^ where F^^ is the electrody- 
namic field tensor and [kp)[LVK\ a tensor that parameter- 
izes Lorentz violation with 19 independent components. 

V.A. Kostelecky and M. Mewes [T? have shown that 
with such an extra term the propagation of light in vac- 
uum can be described in analogy to the propagation of 
light in an anisotropic medium. They have also modeled 
how this anisotropy affects the resonance frequency of a 
linear optical Fabry-Perot resonator. These results can 
be used to model the amplitude coefficients of equations 
^ and ([7| for our experiment. 

The resulting expressions as derived explicitly in |13j 
are given in Table |ll] Ten of the SME parameters are 
linked to birefringence and are restricted to values < 
10~^^ by astrophysical measurements [HI [21]. These 
parameters are assumed to be zero here. Eight of the 
remaining nine SME parameters, grouped into two trace- 
less 3x3 matrices Ke- and Ko+, can then be determined 
from the present measurement: Five are parity even and 
boost independent (Ke_, symmetric) and three are par- 
ity odd and boost dependent (ko-i-, antisymmetric). The 
boost dependent parameters Ko+ lead to an annual phase 
shift of the anisotropy signal due to Earth's orbital revo- 
lution (Table [!!]). Note that since our measurement spans 
more than one year we are indeed able to resolve such an 
annual variation. 

A simultaneous fit of the expressions in Table jllj to the 
2x5 distributions of sidereal modulation amplitudes ob- 
tained from our data, yields estimates on the eight SME 
parameters as summarized in Table jlj Four parameters of 
Ke- and the three boost dependent parameters of /30Ko+ 
(with Earth's orbital boost /J^ = ti^/c = lO""*) feature a 
standard error of - Ix 10^1^, while one parameter, Kg_ , 
shows a slightly increased error bar of 1.7 x 10""'^''. This 
is attributed to the fact that enters the Co compo- 
nent only, which is most prone to the residual systematic 
effects as discussed above in Section fll Al 

All together these limits represent a significant im- 
provement of more than one order of magnitude over the 
results of the best previous experiment by Stanwix et al. 
[1]. They also complement results from a newly emerg- 
ing astrophysical technique, which argues that equivalent 
limits at the 10~^® level can be obtained by analyzing ob- 
servations of ultra-high energy cosmic rays [22) . 
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FIG. 4: Results I. a): Cosine amplitudes C (left) and Sine amplitudes S (right) of a systematic beat frequency modulation at 
2ijJtt for all n — 13384 measurement subsets (10 rotations each), is set to an instant when one of the resonators is oriented 
along the North South axis. b)+c); Amplitudes of a superimposed 2uj(g (11.96 h) (b) respectively cje (23.93 h) (c) modulation 
of C and S respectively, as expected for an anisotropy of c fixed within a sidereal frame. Each point represents the amplitudes 
determined from a 23.93 h set of data. 64 such data sets are included, to is set to an instant when the East- West axis of the 
laboratory coincides with the Y-axis of the adopted SCCEF reference frame. Error bars are omitted for the purpose of clarity 
except for one representative data point, d): Amplitudes Co and So as modeled in equations Q and ([7|, which are most prone 
to constant systematic effects (note the different scale). The mean values and standard errors (shown in red in (b),(c) and (d)) 
of the modulation amplitudes as modeled in equations ([6| and ([7| are: Co = —0.2 ± 11.7, Csi = 15.5 ± 6.0, Cci ~ —9, 9 ± 7.2, 
Cs2 = -1.5±6.6, Cc2 = 4.0±6.6 and So = -10.2±14.4, Ssi = -3.1±6.4, Sci = -2,7±8.1, 5^2 = -5.0±6.2, 5c2 = -5.5±6.7 
(aU values xlO"^*). 
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TABLE I: Estimates on photonic SME parameters obtained 
from this work (one sigma errors). For comparison the limits 
obtained by Stanwix et al. [1] are also given. All values are 
xlO"""^^. /3® = w©/c — 10~* accounts for Earth's orbital 
boost. 
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B. Analysis in the Mansouri-Sexl framework 

We also analyze the data according to the kinematic 
test theory of R. Mansouri and R.U. Sexl j23j, which 
builds on earlier work by H.P Robertson [53]. In this 
test theory a preferred frame is assumed in which the 
speed of light c is isotropic, usually taken to be the cosmic 
microwave background. General, linear transformations, 
using three free parameters a, /?, S, transform from this 
preferred frame to a frame moving at a velocity v. In the 
moving frame an anisotropy of the propagation of light 
then takes the form Ac/c — {(3 + S— |)w^/c^ sin^ 6 where 
6 is the angle between the direction of the propagation 
of light and the direction of v. For a = ^, 13 = ^,6 = 0, 
the generalized transformations reduce to Lorentz trans- 
formations and no anisotropy of c is observed. 

A derivation of the signal amplitudes of equations ^ 
and (It]) in the Mansouri-Sexl framework has been given 
in [13] . The resulting expressions are given in Table III 



Therein we take the velocity of the laboratory relative 
to the CMB as the superposition of the solar system's 
velocity Vc = 370km/s, pointing towards V' = 100° right 
ascension and (/) = —7° declination and the annual mod- 
ulation due to Earth's orbit with — 30km/s. 

Simultaneously fitting these expressions to our data 
yields a value of (/3 + ^ - i) = (4 ± 8) x lO-^^. This is a 
factor of 10 more stringent as compared to the value of 
(9.4 ± 8.1) X 10"" given by Stanwix et al. [Tj. 



III. CONCLUSION 

In conclusion, we have set a limit on an anisotropy of 
the speed of light at a level of Ac/c ~ 1 x 10"-'^'^, which 
allows us to confirm the validity of Lorentz invariance in 
electrodynamics at the 10~^^ level. This accuracy has 
been obtained with optical resonators that feature a rel- 
ative frequency stability of Ai'/vq ^ 1 x 10~^^ in 1 s. The 
final precision could be reached by integrating over more 
than 130 000 rotations relying on a careful suppression of 
systematic effects caused by the turntable rotation. 
Finally, we note that comparable results from a similar 
experiment have been reported after submission of 
this manuscript. 

The relative frequency stability is currently limited by 
thermal noise of the cavity mirrors. Thus, in the longer 
term it should be possible to improve the relative fre- 
quency stability by using cryogenic resonators [111 126] . 
Together with a reasonable improvement in the suppres- 
sion of systematic effects, this would ultimately allow one 
to test for potential violations of Lorentz invariance in 
electrodynamics in the Ac/c ^ 10"^" regime. 
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the Studienstiftung des deutschen Volkes. 
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SME amplitude 

Co: 70 (f^ff - /3(B[(cos77Kf/ + 2 sinrjKf/) cos J^eT + sinfieT]) 
Csi : 71 {-^I- + Ps) [cos 7?Kf/ - sin rjk^^] cos f^ajT) 

Cci : 71 {-af- + [sin ryK^f cos r2®r - ftf/ sin n©r]) 

0,2 ■■ 72 (^fT - I3(B [cos 77^^+^ cos + sin ^^T]) 

Cc2: 72 (i[ftf5 -ftlJ] -/3®[cosr7«f+^cosJ7er-Sr+^sinner]) 



50 : 

5*51 

5*01 



^Ccl 
71 

_ Tar* , 



72 
74 



■^c2 



a 



a2 



TABLE II: Modulation amplitudes according to equations ^ and ^ related to photonic SME parameters. 7o = | sin'^ Xt 
7i = I cos X sin X, 72 = 1(1 + cos^ x)) 73 = —| sinx and 74 = | cosx- Relations are stated to first order in orbital boost. 
/3® = 10~^ is the boost parameter, x ~ 37° is the colatitude of the Berlin laboratory and r) = 23° is the tilt of Earth's axis 
relative to the SCCEF Z-axis. In accordance to the reference frame conventions in [i8| . T = is set to the instant of Earth 
passing vernal equinox. 



RMS amplitude (/? + 5 - |) 

Co'. |7o(— 1 + 3 cos 2(ji) + 2^70 (sin tp cos cos 77 — 2 sin (ji sin 77) cos Q^T + 2^70 cos cos tp sin fi^T 
Csi: —71 sin?/; sin 2(/) — 2^71 (sin <^ cos 77 + sin -(/) cos sin 17) cos ilgT 

Cci: —71 cos 1/; sin 20 — 2^71 cos ip cos sin rj cos fi^T — 2^71 sin (/!>sin Sl^T 

Cb2. —72 sin 2^ cos^ cj> — 2^72 cos t/) cos </> cos 77 cos f7®r — 2^72 sin t/; cos (f> sin fJ^T 

Cc2: —72 cos 2'ip cos^ + 2^72 sin 7/; cos cos 77 cos Q^T — 2^72 cos ip cos sin r2®r 



So: 

SbI 
Scl 
Sb2 
Sc2 





71 

-3a Csi 

71 

?^a2 



TABLE III: Modulation amplitudes according to equations ^ and Q related to the RMS parameter (/3 + <5 — 5). 70 = 

ijsin2Xi7i — 5 sinx cosx, 72 = i(l + cos'^x):73 = ~|sinx and 74 = |cosx- X = 37° denotes the laboratory colatitude, 
r) = 23° the tilt of Earth's axis relative to the SCCEF Z-axis. tp — 100° right ascension and <j) — —7° declination denote the 
direction of the solar systems velocity relative to the cosmic microwave background. Terms varying with Earth 's orbital motion 
are suppressed by v^/vc ~ 0.08. Aiso here T = is set to the instant of Earth passing vernal equinox. 
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